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Abstract 

We review recent progress in the study of transport properties of inter- 
acting electrons subject to a disordered potential which is strong enough 
to localize all single-particle states. This review may also serve as a guide 
to the recent paper by the authors [Annals of Physics (2006), in press]. 
Here we skip most of the technical details and make an attempt to discuss 
the physical grounds of the final-temperature metal-insulator transition 
described in the above-mentioned paper. 

1 Introduction 

Transport properties of conducting materials at low temperature T are deter- 
mined by an interplay between the interaction of the itinerant electrons with 
each other and the quenched disorder which creates a random potential act- 
ing on these electrons. In the absence of the electron-electron interaction the 
most dramatic phenomenon is Anderson localization (Anderson, 1958) - the dc 
electrical conductivity a can be qualitatively different depending on whether one- 
particle wave functions of the electrons are localized or not. In the latter case 
er(T) has a finite zero-temperature limit, while in the former case a(T) vanishes 
when T — > 0. Therefore, Anderson localization of electronic states leads to the 
Metal to Insulator Transition at zero temperature. 

When discussing zero temperature conductivity tr(0), we need to consider 
only electronic states close to the Fermi level. The conductivity becomes finite 
at any finite temperature provided that extended states exist somewhere near 
the Fermi level. It is commonly accepted now that localized and extended states 
in a random potential can not be mixed in the one-electron spectrum and thus 
this spectrum in a general case is a combination of bands of extended states 
and bands of localized states. A border between a localized and an extended 
band is called mobility edge. If the Fermi level is located inside a localized band 
and inelastic scattering of the electrons is completely absent, the conductivity 
should follow Arrhenius law <j(T) oc exp(— E c /T), where E c is the distance from 
the Fermi level to the closest mobility edge. Another common belief following 
from the scaling theory of Anderson localization (Thouless, 1977; Abrahams et 
al, 1979), is that in low dimensionality d, namely at d = 1,2 all states are 
localized in an arbitrarily small disorder, while for free electrons (no periodic 
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potential) E c > is finite at d = 3. 1 It means that without inelastic processes 
Cd=i,2(T) — 0, while for <r d=3 (T) one should expect the Arrhenius law. Note that 
for electrons in a crystal within a given conduction band the latter conclusion is 
not always correct - strong enough disorder can localize the whole band. 

As soon as inelastic processes are included, the situation becomes more 
complicated. In particular, electron-phonon interaction leads to the mechanism 
of conductivity known as hopping conductivity (Fritzche, 1955; Mott, 1968a; 
Shklovskii and Efros, 1984) - with an assistance of phonons, electrons hop be- 
tween the localized states without being activated above the mobility edge. As 
a result, <j(T) turns out to be finite (although small) at arbitrarily low T even 
when all one-electron states are localized. 

Can interaction between electrons play the same role and cause the hopping 
conductivity? This question was discussed in literature for a long time (Fleishman 
and Anderson, 1980; Shahbazyan and Raikh, 1996; Kozub, Baranovskii, and 
Shlimak, 2000; Nattermann, Giamarchi, and Le Doussal, 2003; Gornyi, Mirlin, 
and Polyakov, 2004) and no definite conclusion was achieved. The problem is that 
although the electric noise exists inside the material with a finite ac conductivity 2 
the "photons" in contrast with phonons become localized together with electrons. 

In a recent work (Basko, Aleiner, and Altshuler, 2006) we have demonstrated 
that electron-electron interaction alone cannot cause finite conductivity even 
when temperature is finite, but small enough. In the absence of phonons and 
extended one-electron states conductivity of a system of interacting electrons 
vanishes exactly below some critical temperature T c . At the same time, at high 
temperatures T > T c the conductivity a(T) is finite. It means that at T = T c the 
system of interacting electrons subject to a random potential undergoes a genuine 
phase transition that manifests itself by the emerging of a finite conductivity! 

This transition can be thought of as many-body localization - it applies to 
many-body eigenstates of the whole system. This localization occurs not in the 
real space, but rather in the Fock space. This fact does not affect the validity 
of the concept of mobility edge. In fact, the existence of the "metallic" state at 
T > T c implies that the many-body states with energies £ above £ c are extended. 
One can estimate the difference between £ c and the energy of the many-body 
ground state £q as £ c — £o ~ TJ\F(T), where J\T{T) is the total number of one- 
particle states in the energy strip of the width T. Note that the existence of 
the extended many-body states above the mobility edge does not contradict the 
fact that below T c there is no conductivity - in contrast with the case of one- 
particle localization there is no Arrhenius regime since £ c — £q turns out to be 
proportional to the volume of the system, i.e., is macroscopically large (see Sec.0] 
below for more details). 

1 For d = 1 this statement was proved rigorously both for one-channel (Gertsenshtein and 
Vasil'ev, 1959; Berezinskii, 1973) and multi-channel (Efetov and Larkin, 1983; Dorokhov, 1983) 
disordered wires. 

2 In this paper we mostly focus on dc conductivity. As to ac conductivity, it never vanishes, 
because at any frequency density of resonant pairs of states is finite. 
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In order to avoid possible misunderstanding we would like to emphasize that 
we focus only on the inelastic collisions between the electrons, i.e., on creation 
or annihilation of real electron-hole pairs. There are other effects of electron- 
electron interactions which can be understood as renormalization of the one- 
particle random potential by the interaction. Being temperature-dependent, this 
renormalization leads to a number of interesting effects, such as the interac- 
tion corrections to the density of states and conductivity in disordered metals 
(Altshuler and Aronov, 1985). On the insulating side of the one-particle localiza- 
tion transition similar effects cause the well-known Coulomb gap (Shklovskii and 
Efros, 1984) which reduces hopping conductivity. On the other hand, this is just 
a correction to the time-independent random potential. As such, it can maybe 
shift the position of the many-body Metal to Insulator transition, i.e., renormal- 
ize T c , but is unable to destabilize the insulating or metallic phases. From now 
on we will simply neglect all elastic (Hartree-Fock) effects and concentrate on 
the real inelastic electron-electron collisions. 

Localization of the many-body states in the Fock space has been discussed 
by Altshuler et al. (1997) for the case of zero-dimensional systems with finite, 
although large, number of electrons. In this paper the authors proposed an ap- 
proximate mapping of the Hamiltonian of a metallic grain with large Thouless 
conductance g and moderate interaction between the electrons to the one-particle 
Hamiltonian on a lattice with the topology of the Cayley tree and an on-site dis- 
order. The latter problem has an exact solution (Abou-Chacra, Anderson, and 
Thouless, 1973; Efetov, 1987) that exhibits the localization transition. In terms 
of interaction electrons this transition means that one-particle excitation states 
below certain energy are quite close to some exact many-body excitations. As to 
the one-particle excitations with energies higher than the critical one, its wave 
function can be viewed as a linear combination of a large number of the many- 
body eigenstates. 

For an infinite system (d > 0) the situation is more complex, and Cayley 
tree approximation is hard to justify. Nevertheless, a consistent analysis of a 
model with weak and short range interaction to all orders of perturbation theory 
enabled us to analyze the many-body localization transition and to demonstrate 
that both the metallic state at high temperatures and the insulating state at low 
temperatures are stable and survive all higher loop corrections to the locator 
expansion. Therefore, the existence of the transition is proved on the physical 
level of rigor. 

It should be noted that such an insulating state that is characterized by ex- 
actly zero conductivity is quite different from all other known types of insulators. 
For example, Mott insulator is believed to have finite, though exponentially small 
conductivity at finite temperatures. 

The present text represents a shortened version of the paper by Basko, 
Alciner, and Altshuler (2006), hereafter referred to as BAA paper. We omit 
most of the technical details (for which the reader will be referred to specific 
sections of the BAA paper), and stress the key ideas. 
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The remainder of the paper is organized as follows. In Sec.|21we briefly review 
some well-known facts about electric conduction in Anderson insulators and 
pose the problem. Sec. El represents a sketch of the solution whose details are 
given in BAA paper. We discuss the model for interacting localized electrons in 
Sec. 13.11 and the corresponding Fock space picture in Sec. 13.21 In Sec. 13.31 we 
show the formal way to characterize metallic and insulating phases. In Sec. 13.41 
we introduce the main approximation used in the calculation (self-consistent 
Born approximation), and discuss its validity. The existence of the metallic state 
at high temperatures and its properties are discussed in Sec. 13.51 Sec. 13.61 is 
dedicated to the proof of existence of the insulating phase at low temperatures; 
the value of the transition temperature is obtained as the limit of stability of 
the insulating phase. In Sec. 21 we discuss the macroscopic implications of the 
problem, introducing the concepts of many-body localization and many-body 
mobility edge. Finally, in Sec. we summarize the results and present an outlook 
of the future developments. 

2 Background and formulation of the problem 

2.1 Non-interacting electrons in disorder potential 

Let us briefly review the basic concepts developed for the problem of one-electron 
wave functions in a disordered potential in d dimensions. Depending on the 
strength of the disorder potential, a wave function </> a (r) of an eigenstate a with 
the energy £ Q can be either localized or extended: 



Here Q oc is the localization length which depends on the eigenenergy £ Q , and 
Q is the volume of the system. Each localized state is characterized by a point 
in space, p a , where l^a^l 2 reaches its maximum, and an exponentially falling 
envelope. Extended states spread more or less uniformly over the whole volume 
of the system. Localized and extended states cannot coexist at the same energy, 
and the spectrum splits into bands of localized and extended states. The energies 
separating such bands are known as mobility edges. For free electrons in d > 3 
disorder potential leads to only one mobility edge £i , so that 



If a finite mobility edge J2J) exists and the Fermi level cf lies in the band of local- 
ized states, the conductivity is determined by the exponentially small occupation 
number of the delocalized states 



|^(r)| 2 cx 




, localized; 
extended. 



(1) 



£ a < £\ : localized; 
£ Q > E\ : extended. 
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In this paper we are interested in transport properties of the systems where 
all single-particle states are localized, and thus without many-body effects a = 
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at any temperature. It is well established now that the mobility edge usually does 
not exist for one- and two-dimensional systems, and all single-particle states are 
indeed localized for an arbitrarily weak disorder. Such a situation can arise for 
a large d as well, if the bandwidth is finite and disorder is sufficiently strong. 

2.2 Role of inelastic processes and phonon- assisted hopping 

As long as all single-particle states are localized, transport occurs only because 
of inelastic processes, which transfer electrons between different localized eigen- 
states. At this stage we introduce the main energy scale of the problem: the typ- 
ical energy spacing between states whose spatial separation does not exceed C,i oc , 
so that there is overlap between their wave functions: 

*< = ^' (4) 

Woe 

where v is the one-particle density of states per unit volume. 

The conductivity is, roughly speaking, proportional to the rate of the transi- 
tions between different localized states, which can be called inelastic relaxation 
rate T. Obviously, at T = inelastic processes disappear, so regardless of the 
mechanism of inelastic relaxation the conductivity must vanish: 

lira a(T) = . (5) 

The question is how <r(T) approaches zero for each particular mechanism. 

When phonons are the main source of inelastic scattering, the answer is given 
by Mott's variable range hopping formula (Mott, 1968a) 3 



cr(T) =o- [ — \ exp 



(6) 



where Co and a are constants. According to Eq. (JBJ, er(T) remains finite as long 
as T ^ 0. The reason for cr(T) ^ is that for any pair of localized states one 
can always find a phonon whose frequency exactly corresponds to their energy 
mismatch, at low temperature one should just wait long enough. 

The same type of <r(T)-dependence would result from the coupling of elec- 
trons with any delocalized thermal bath whose energy spectrum is continuous 
down to zero energy. The specific nature of the bath at most affects the power- 
law prefactor. On the contrary, the stretched exponential factor is universal; it 
originates from the counting of electronic states available for the transition, and 
does not depend on the specific scattering mechanism. 

2.3 Inelastic relaxation due to electron- electron interaction 

Now let us assume that there is no external bath coupled to electrons, but some 
electron-electron interaction is present. What will be the dependence of cr{T)l 

3 Here we do not consider effects of the Coulomb interaction which are known to modify the 
power of temperature in the exponent (Shklovskii and Efros, 1984). 
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In line with the discussion of the previous subsection, the question should be 
posed as follows: do electron-hole pairs themselves provide a suitable bath in a 
localized system, thus validating Eq. JBJ? 

One possible answer is "yes" . Indeed, recall Mott formula for the low-tempera- 
ture dissipative ac conductivity a(u>) in a localized system (Mott, 1968b): 

(j(oj) = <7i ln t^t ■ (7) 



According to fluctuation-dissipation theorem, at finite temperature electromag- 
netic fluctuations of a finite spectral density should be present, and they might 
serve as a bath. The problem with this argument is that Eq. Q is the spatial 
average of the conductivity over the whole (infinite) volume. For each given re- 
alization of disorder, excitations determining <t(oj) from Eq. Q are localized, 
and although the total volume is infinite, the spectrum of electron-hole pairs is 
effectively discrete. 

The crucial point is that as long as electron-electron interaction is local in 
space (here we do not consider long-range interactions), it effectively couples 
electronic states only within the same localization volume, where the spectrum 
of electronic states is effectively discrete. The following sections are dedicated to a 
systematic discussion of this problem which was first pointed out by Fleishman 
and Anderson (1980). The conclusion can be stated as follows: electron-hole 
excitations can cause finite conductivity only if the temperature of the system 
exceeds some critical value. At lower temperatures a(T) vanishes exactly. 

3 Finite-temperature metal-insulator transition 

3.1 Matrix elements of electron- electron interaction between localized states: 
essential features of the model 

For simplicity we consider a system of spinless electrons and assume that electron- 
electron interaction is weak and short-range: 4 

V(n - rfe) = £ 6(n - f 2 ) , (8) 

where A <C 1 is the dimensionless interaction constant, v is the one-particle 
density of states per unit volume. 

In the basis of localized single-particle eigenstates the hamiltonian corre- 
sponding to the pair interaction potential JSJl takes the form 



H = £ a + \ V a (3 7 s clcl&ycs (9) 



4 Interaction proportional to S(f\ — Fa) in the strict sense is equivalent to no interaction for 
spinless electrons, considered here, due to the Pauli principle. Here, by writing &{r\ — r^) we 
only mean that the range is much smaller than the electron mean free path, so it is not a true 
5-function. 
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Consider the structure of the matrix elements Va^s- Since V(f) is short-range, 
they decrease exponentially when the spatial separation between the states in- 
creases, the characteristic scale being the localization length Q oc . In addition to 
this spatial suppression, the matrix elements decrease rapidly when the energy 
difference, say £ Q — £ 7 , increases exceeding the level spacing 5^. This occurs be- 
cause the localized wave functions oscillate randomly, and the bigger the energy 
difference, the weaker are these random oscillations correlated (Altshulcr and 
Aronov, 1985). Provided that the restrictions 

\r a - rp\ < doc, \r a - f y \ < Cioc, \rp - r 7 | < Q oc , etc., (10) 
|&-&|,|fr-erl<<fc or |£«-ey|,|&-6|<5 C . ( n ) 

are fulfilled, we have |V Q( 3 7< 5| ~ XS^. 

There are several ways to model these essential properties of V a p 7 s- In BAA pa- 
per a specific model was adopted; essentially, the space and energy dependences 
of the matrix elements were replaced by simple rectangular cutoffs (see Sec. 3 of 
BAA paper for details). 

3.2 Many-electron transitions and Fock space 

We wish to note that all the discussion of this subsection is not conceptually 
new. In fact, it is just a generalization of the arguments of Altshulcr et al. (1997) 
to an infinite system. 

Conventionally, an elementary inelastic process is a decay of one single- 
particle excitation (an electron occupying a state a) into three single-particle 
excitations - a hole in the state /3 and two electrons in the states 7 and 6. Such 
a decay can be described differently: one can say that the hamiltonian couples 
the single-particle excitation with the three-particle excitation by the matrix ele- 
ment Va/3-ys- Further action of the interaction hamiltonian produces five- particle 
excitations, seven-particle excitations, etc.: 

£a -> £y+&-& - . (12) 

If on each stage the coupling is strong enough (i. e., the matrix element is of 
the same order or larger than the corresponding energy mismatch), the single- 
particle excitation indeed decays irreversibly into all possible many-body states. 
In other words, exact many-body eigenstates become delocalized in the Fock 
space. If, oppositely, three-particle states contribute only a weak perturbative 
admixture to the one-particle state, the contribution from five-particle states is 
even weaker, etc., the initial electron will never decay completely. One can say 
that it is localized in the Fock space. 

Another way to visualize the inelastic relaxation is to look at the energy 
structure of the quasiparticle spectral function: 

^a(e) = Yl K*fcl4l*o>| 2 S(e + E Q - E k ). (13) 
k 
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A A (1) (e) 



(a) 



A A (3) (e) 



(b) 




A A f7) (8) 



(d) 



IiLll 



^2 ^3 '4 '5 



E +1 +1 +t -t -E -t e 

'2 ?3 > ^5 '6 ^7 



Fig. 1 . A schematic view of the contributions to the perturbative expansion (|14fl 
of the spectral function A a (e) from one-particle, three-particle, five-particle, 
and seven-particle excitations. 

Here and "Jfc are many-body eigenstates (^o is not necessarily a ground 
state), E and E k are the corresponding energies. Basically, A a (e) shows how 
the single-particle excitation on top of a given eigenstate is spread over other 
many-body eigenstates of the system. 

One can represent the spectral function in the form of expansion in powers 
of the interaction constant A (Fig.^): 



oo 

A a (e) = yX 2n A^(e). 



n=0 



(14) 



Aa\e) corresponds to the bare quasiparticle peak, 5(e — £ Q ) . Linear in A term 
represents a random Hartree-Fock shift of the energy £ Q which is already assumed 
to be random. This linear term is of no interest to us, so it is not included in the 
expansion l|14l) . The A 2 term corresponds to the contribution of three-particle 
excitations. The number of these excitations is effectively finite due to the re- 
strictions HI U| l. TO ; so A( 3 \e) is a collection of <5-peaks at three-particle energies 
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£7 + £5 — C/3- Again, A terms correspond only to peak shifts and are omitted. 
The five-particle contribution A^ 5 \e) is again a collection of <5-peaks, however, 
spaced more closely than three-particle peaks. This represents the general rule: 
the more particles are involved, the higher is the density of the peaks. At the 
same time, contributions of many-particle processes are suppressed due to the 
smallness of A. What is crucial to us, is the result of this competition at n — > 00: 
will the empty spaces between peaks be filled making A a (e) a continuous func- 
tion, or the process will be suppressed by the powers of A and A a (e) will remain 
a collection of 5-peaks? 

Thus, the problem of inelastic relaxation of a quasiparticle due to electron- 
electron interaction is related to the problem of localization or derealization 
of excitations in the many-electron Fock space. The simplest model describing 
localization-delocalization physics is Anderson model. Anderson (1958) consid- 
ered a tight-binding model on a d-dimensional lattice. The coupling between sites 
is nearest-neighbor only with a fixed matrix element V. The on-site energies are 
assumed to be random and uncorrelated, with some typical value denoted by W. 

It is thus tempting to identify the many-electron hamiltonian with the 
Anderson hamiltonian on a certain lattice, whose sites correspond to many- 
particle excitations (QJ). The approximate rules of correspondence then should 
be the following: 

• V — > A<5^ - typical value of the coupling matrix element; 

• W — > 8( - typical energy mismatch in each consecutive virtual transition, 

If a +& ~ fy - &| ~ % 

• finally, the coordination number 2d — > T/5q. This represents the number 
of three-particle excitations to which a given single-particle excitation is 
coupled, with the energy mismatch not exceeding 5q. The value T/5q is 
obtained as the product of the number of electrons (3 within the same 
localization volume, available for the collision with the probe electron a 
(~ T/8^), and the number of ways to distribute energy allowed by the 
restriction (jlip. which is ~ 1. 

According to Anderson (1958), localization-delocalization transition occurs at 

Vd W XT 1 

^ ln V™ A (15) 

For the interacting localized electrons this would imply that below some tem- 
perature T c the inelastic relaxation is frozen, so that the conductivity vanishes 
exactly; above T c some inelastic relaxation is taking place , and the conductivity 
is finite. This corresponds to a finite-temperature metal-insulator transition. 

This analogy, however, should be used with caution. Anderson's result, in 
fact, can be sensitive to the structure of the lattice; possible sources of problems 
are listed in Sec. 2.2 of BAA paper. One can develop a systematic approach 
to the problem, based on the diagrammatic technique for interacting electrons 
(Basko, Aleiner, and Altshuler, 2006). Below we headline the basic ideas of this 
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Fig. 2. (a) Schematic energy dependence of the quasiparticle decay rate T(e) 
in the metallic (dashed line) and insulating (solid line) phases for a given 
realization of single-particle levels {£ Q }. (b) The corresponding distribution 
functions P(T) for T(e) at a given energy e in the metallic (dashed line) and 
insulating (solid line) phases. 

approach to justify the assumption that Eq. I|15|) indeed correctly determines the 
point of the metal-insulator transition. 

3.3 Statistics of the transition rates 

We have to admit, that all the discussion of this subsection is not conceptually 
new either. It rather represents a generalization of the arguments of Anderson 
(1958) to a many-body problem. 

As discussed in previous sections, the focus of the problem is the inelastic 
quasiparticle relaxation, which is represented by the imaginary part of the single- 
particle self-energy: 

r Q ( e )=Im^(e). (16) 

We stress that r a (e) is a random quantity, as it depends on the positions of all 
single-particle levels {£,3} and all occupation numbers {np}. One has no other 
choice but to perform statistical analysis of this random quantity. 

How to distinguish between metallic and insulating phases within a statistical 
framework? It is clear that positions of the peaks in r(e) for the insulating regime 
wander randomly with the variation of random energies £ a from Eq. JJjJ . Due to 
this variation, the ensemble average of the decay rate (r(e)) is the same in both 
phases and can not be used for the distinction. In fact, one has to investigate 
the whole distribution function P(T). 

To understand the difference in the behavior of P(T) in the two phases, it is 
instructive to start with the behavior of r a (e) for a given realization of disorder, 
Fig. E^,. Deep in the metallic phase r Q (e) is a smooth function of energy, so its 
distribution P(T) at a given energy is a narrow gaussian. In the insulating phase 
r Q (e) is given by a sequence of infinitely narrow ^-peaks, so an arbitrarily chosen 



Fimte-temperature metal-insulator transition 



11 




a 




a 



+ 




i 



Fig. 3. Diagrammatic representation of self-consistent Born approximation. 

value of e falls between the peaks with the probability 1, giving T = 0. If e hits 
a (5-peak, then r = oo, which happens, however, with zero probability. 

One can deal with this uncertainty by introducing an infinitesimal imaginary 
energy shift (damping) t], and consider T a (e + irj). Physically, it may be viewed 
as an infinitesimally weak coupling to a dissipative bath (e. g., phonons). Clearly, 
in the metallic phase it has no effect, while in the insulating phase this damping 
broadens (5-peaks into lorentzians of the width 77, which leads to appearance of 
the tails. Now, even if the energy e falls between the peaks, V will have a finite 
value proportional to rj. As a result, the distribution function P(T) will have the 
form sketched in Fig. [2t>- Now, having calculated the distribution function for a 
small but finite 77, one can distinguish between the metallic and the insulating 
phases according to 



Here fl is the total volume of the system. We emphasize that the order of limits 
in Eq. Ijl7(l cannot be interchanged, since in a finite closed system the spectrum 
always consists of discrete <5-peaks. 

3.4 Self- consistent Born approximation 

Our main object of interest is the probability distribution function of the decay 
rate P{T). Its calculation is performed in two stages: (i) we find T a (e) for a given 
realization of disorder, (ii) we calculate its statistics. This subsection is dedicated 
to the first task. 

We intend to describe both metallic and insulating regimes. In the latter 
regime relaxation dynamics is absent, the system never reaches thermal equilib- 
rium, and temperature itself is ill-defined. Therefore, the only appropriate formal 
framework is the non-equilibrium formalism of Keldysh (1964) (see Sec. 4.1 of 
BAA paper for details). 




(17) 



12 



D. M. Basko, I. L. Aleiner, and B. L. Altshuler 



Our approach for calculation of r a (e) is the self-consistent Born approxima- 
tion (SCBA), shown diagrammatically in Fig- El an d given by 



Here n a — 0, 1 is the fermion occupation number of the single-particle state a. 
It is crucial to realize that one cannot replace n a with its average equilibrium 
value, corresponding to Fermi-Dirac distribution. The basis in the Fock space is 
formed by Slater determinants corresponding to n a = 0, 1 only, and T represents 
the rate of a transition between two such basis states with different sets of {n a }. 
Temperature enters through total energy of the state, which is determined by 
all {n a } and is proportional to the volume of the system (this issue will be 
discussed in more detail in Sec. El of the present paper). 

Iterations of SCBA generate self-energy diagrams which have one common 
property: they describe decay processes where the number of particles in the 
final state is maximized for each given order in A, thus maximizing the phase 
space available for the decay. Vice versa, each diagram satisfying this condition 
is taken into account by SCBA. Let us briefly discuss contributions which are 
neglected by Eqs. (|18J) . (|19|1 : for a detailed discussion see Sec. 7 of BAA paper. 

• Eqs. I|18l) . (|19fl completely ignore the real part of the self-energy, ReS. In 
most of the terms effect of Re £ reduces to random uncorrelated corrections 
to already random energies, and appear to be completely negligible. In a 
more accurate approximation, however, some statistical correlations are 
present. They renormalize the numerical prefactor in the expression for 
the critical temperature, see Sec. 7.3 of BAA paper. 

• Generally, quantum mechanical probability of a transition is given by a 
square of the total amplitude, the latter being a sum of partial amplitudes. 
In fact, Eqs. (|18fl . (|19fl correspond to replacing the square of the sum by 
the sum of the squares. Approximation which neglects the interference 
terms can be justified in the same way as in the Anderson model of a high 
dimensionality d e ff ~ T/Sq ~ 1/A ^> 1, see Sec. 7.2 of BAA paper. 

• Finally, diagrams generated by SCBA correspond to taking all n-particle 
vertex function in the leading approximation in A. One can show that for 
A <C 1 vertex corrections are small, see Sec. 7.1 of BAA paper. For A ~ 1 
one would have to introduce the full n-particle vertex analogously to how 
it is done in conventional Fermi liquid theory for n = 2 (Landau, 1958; 
Eliashberg, 1962). In this case the result for the transition temperature will 
be determined not by the bare interaction constant, but by the statistics of 
the full vertex functions. Existence and regularity of these vertex functions 




0,i,s 



LO) X 



x [np(l - n 7 )(l - n s ) + (1 - n /3 )n 7 n (5 ] , 

1 r a (e) 

* (e-6*) 2 + r2( £ )- 



(18) 



(19) 
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are basically equivalent to the assumption that in the absence of disorder 
the interacting system is a Fermi liquid. 

3.5 Metallic phase 

Starting from the self-consistent equations JTSJ, (|19|l . one can straightforwardly 
verify (see Sec. 5.1 of BAA paper) that P(F) can be well approximated by 
gaussian distribution with the average and dispersion 

(r) ~ a 2 t, (r 2 ) - (r) 2 ~ a 2 <5 2 , (20) 

provided that 

V(r 2 ) _ (fp < ( F ) ^ r>T (»)_^. ( 2 i) 

A 

According to the arguments of Sec. 13.31 this is characteristic of the metallic 
phase. One should not think, however, that this automatically means that the 
system has the same transport properties as conventional metals (which would 
mean that the conductivity is given by the Drude formula) . The system conducts 
in the Drude regime only when the inelastic processes completely destroy the 
localization, the latter manifesting itself as the weak localization correction to 
conductivity (Altshuler and Aronov, 1985). This occurs when the discrete levels 
are completely smeared: 

(r)»<5 c O T»T( e ')~^. (22) 

It turns out that for A ^ 1 there is a parametric range of temperatures: 

^<T<^|. (23) 

In this interval electron-electron interaction is already sufficient to cause inelastic 
relaxation, however, the localized nature of the single-particle wave functions 
remains important. Conduction in this regime can be viewed as an analog of 
hopping conduction discussed in Sec. 12.21 in the sense that electrons themselves 
indeed provide a good bath. However, as T ^> 5^, there is no exponential factor 
in the temperature dependence cr(T). This dependence can be obtained from the 
kinetic equation describing electron transitions between localized states, and is 
given by a power law, a(T) oc T a , where a can be model-dependent. Also, in this 
temperature range Wiedemann-Frantz law can be violated. For further details 
the reader is referred to Sec. 5.2 of BAA paper. We only note that this regime 
is somewhat analogous to the phonon-assisted conduction discussed by Gogolin, 
Mel'nikov, and Rashba (1975). 

3.6 Insulating phase 

Self-consistent equations 118|) , 119|) represent a system of nonlinear integral equa- 
tions whose coefficients are random due to randomness of level energies {£, a } and 
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occupation numbers {n a }. Apparently, for r\ — these equations have a solu- 
tion T a (e) = 0, corresponding to the insulating phase. One must check, however, 
whether this solution is stable with respect to an infinitesimal damping ?y. 

To perform the standard linear stability analysis of Eqs. 118|l . I|19|) we linearize 
Eq. (Uni as 

A a (e) = S(e - + i + 0(T 2 ) , (24) 

in complete analogy with Abou-Chacra, Anderson, and Thouless (1973), substi- 
tute this into Eq. Q18)). and obtain a linear integral equation: 

V (A „+X-\V ,2 2r 7 (e + &-&) + r>(g 7 +&-e) 



x [^(l - n 7 )(l - 715) + (1 - np)n~jn$] . (25) 



The solution of this equation may be sought in the form of a perturbations series 

oo 

r Q (e) = £rM(e), (26) 

where is of the order \V\ 2n ~ (\S() 2n and is obtained after n iterations of 
Eq. H25f) starting from r' ^ = rj. Each term in this expansion can be calculated 
explicitly and its statistics can be determined (details of this rather cumbersome 
calculation are given in Sec. 6 of BAA paper). The resulting probability dis- 
tribution function is controlled by a single parameter 7„, which determines the 
typical scale of the random quantity r^: 

P„(r(")) = e --/ r( " , ./-^3, 7 „ = C^A 2 ", A = C 2 ^lni (27) 

y [r(»)] k x 

Here C\ ~ 1 and C2 ~ 1 are model-depndent numerical constants (Eq. (172) of 
BAA paper contains their values for the specific model adopted there). 

The behavior of this distribution function at n — + 00 is qualitatively different, 
depending on whether A in Eq. 1)2 7JI is smaller or larger than 1, as illustrated by 
Fig-E] In the first case the typical scale 7„ — > as n — -> 00, and the distribution 
function for the total Y is always concentrated around r ~ r\ which tends to zero 
as r\ — » 0. This perfectly matches the insulating behavior described in Sec. 13.31 
and shown in Fig. 03. 

Oppositely, for A > 1 the scale j n — > 00; we emphasize that the limit n — > 00 
should be taken prior to 77 — > 0. This means that the distribution function 
P(T) does not shrink to S(T) as rj — > 0; in fact, it cannot be determined from 
the linearized equation l|25l) . the full self-consistent problem (|18fl . (|19|) should 
be solved. The divergent linear solution signals the instability of the insulating 
state and the onset of the metallic state. 
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Fig. 4. Plot of the distribution functions P n (T^) as given by Eq. I|27|) for 
n = 1, 2, 3 for T/T c = 0.7 (solid lines) and T/T c = 1.3 (dashed lines). 

These arguments enable us to identify the temperature at which A = 1 with 
the temperature of the metal-insulator transition: 



4 Metal-insulator transition and many-body mobility edge 

In the previous section we considered the decay of a quasiparticle excitation, and 
found that possibility or impossibility of such decay may be viewed as dereal- 
ization or localization of this excitation in the many-electron Fock space. In this 
section we discuss macroscopic implications of this picture. How can the very 
notion of localization be applied to many-body states? 

Consider a many-body eigenstate \^k) of the interacting system, with the 
corresponding cigenenergy E^. In the coordinate representation, the many-body 
wave function "ffc {{fj^—^ depends on the coordinates of all N particles in the 

system. Let us create an electron-hole pair on top of \^k)- The resulting state, 
which is not an eigenstate of the system, can be expanded in terms of other 
eigenstates: 



It is possible that the number of terms contributing to the sum is effectively 
finite, i.e. 



C 2 Aln(l/A) ' 



(28) 




(29) 
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-l 

< oo. (30) 

This corresponds to insulating or localized many-body state; excitation can not 
propagate over all states allowed by the energy conservation. 

The opposite case, when expansion (|29[) contains an infinite number of eigen- 
states 

= oo, (31) 

corresponds to metallic or extended many-body state. 

Developed metallic state is formed when expansion l|29(l involves all the eigen- 
states with close enough energies: 

|C^'| 2 oc «6(E k +u aP -E k ,y>, (32) 

where 5-function should be understood in the thermodynamic sense: its width, 
although sufficiently large to include many states, vanishes in the limit £1 — > 
oo. Only in this regime, which may also be called ergodic many-body state, the 
electron-electron interaction can bring the system from the initial Hartree-Fock 
state to the equilibrium corresponding to spanning all the states permitted by 
the energy conservation. In this case, the averaging over the exact many-body 
eigenfunction is equivalent to averaging over the microcanonical distribution, 
and temperature T can be defined as a usual Lagrange multiplier. It is related 
to E k by the thermodynamic relation: 

T 

E k -E = JcvtTjdTu (33) 
o 

where Eo is the ground state energy, and CV(T) oc f2 is the specific heat. 

The temperature of the metal-insulator transition, found above, in fact, de- 
termines the extensive many-body mobility edge £ c oc il. In other words, (i) states 
with energies E k — Eq > £ c are extended, inelastic relaxation is possible, and 
the conductivity o~ k = o~(E k ) in this state is finite; (ii) states with energies 
E k — Eq < £ c are localized and the conductivity a(E k ) = 0. 

Let us now assume that the equilibrium occupation is given by the Gibbs 
distribution. One could think that it would still imply the Arrhenius law @ 
for the conductivity. However, this is not the case for the many-body mobility 
threshold. In fact, in the limit fl — > oo 

a(T) = 0; T < T c , (34) 

where the critical temperature is determined by Eq. I|33|l : 



lim 

V^oo 



E 

. k' 



c, 



kk' 

a0 



lim 



E 

. k' 



Ckk' 

I a/3 



/ ' dT 1 C v (T 1 )=£ c . 
Jo 



(35) 
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Therefore, the temperature dependence of the dissipative coefficient in the system 
shows the singularity typical for a phase transition. 

To prove the relations (|34[) . (|35[) we use the Gibbs distribution and find 

^ / °° dEe s ^- E / T a{E) 

<*{T) = 2-.Pk0{E k ) = J°° dEeS{E )-E/T ' 

where the entropy S(E) is proportional to volume, and E is counted from the 
ground state. The integral is calculated in the saddle point or in the steepest 
decent approximations, exact for Q — > oo. The saddle point E(T) is given by 



dS_ 

dE 



1 

f 



E=E(T) 

Taking into account a(E) = for E < £ c we find 
a(T)=a[E(T)], E(T) > £ c ; 
<r(T) ocexpf- g^S V E(T) < £ c 

As both energies entering the exponential are extensive, E(T), £ c oc f2, we obtain 

To be able to establish the thermal equilibrium in such insulating state the 
system should be coupled to an external bath (i.e., phonons). The presence of the 
finite electron-phonon interaction (as phonons are usually delocalized), smears 
out the transition, and <r(T) becomes finite for any temperature. Nevertheless, if 
electron-phonon interaction is weak, the phenomenon of the many-body metal- 
insulator transition manifests itself as a sharp crossover from phonon induced 
hopping at T < T c to the conductivity independent of the electron-phonon cou- 
pling at T > T c . 



5 Conclusions and perspectives 

We have considered the inelastic relaxation and transport at low temperatures 
in disordered conductors where all single-particle states are localized, and no 
coupling to phonons or any other thermal bath is present. The main question 
is whether electron-electron interaction alone is sufficient to cause transitions 
between the localized states producing thereby a finite conductivity. The an- 
swer to this question turns out to be determined by the Anderson localization- 
delocalization physics in the many-particle Fock space, and is summarized on 
Fig. El 

It should be emphasized that the many-body localization, which we discuss in 
this paper, is qualitatively different from conventional finite temperature Metal 
to Insulator transitions, such as formation of a band insulator due to the struc- 
tural phase transition or Mott-Hubbard transition. In these two cases, at a cer- 
tain temperature T* a gap appears in the spectrum of charge excitation (Mott 
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Ao-(T) k k k 




Fig. 5. Schematic temperature dependence of the dc conductivity <j(T). Below 
the point of the many-body metal-insulator transition, T < T c , no inelastic 
relaxation occurs and a(T) = 0. Temperature interval T > T {in ^ > T c 
corresponds to the developed metallic phase, where Eq. Q21JI is valid. At 
T 3> the high-temperature metallic perturbation theory (Altshuler and 
Aronov, 1985) is valid, and the conductivity is given by the Drude formula. 

insulator) or all excitations (band insulator). In both cases the conductivity re- 
mains finite although exponentially small as long as T > 0. This is not the 
case for many-body localization, which implies exactly zero conductivity in the 
low-temperature phase. 

Is the many-body localization a true thermodynamical phase transition with 
corresponding singularities in all equilibrium properties? This question definitely 
requires additional studies, however, some speculations can be put forward. 
The physics described in the present paper is associated with the change of 
the characteristics of the many-body wavefunctions. It is well known that for 
non-interacting systems localization-delocalization transition does not affect the 
average density of states, i.e., it does not manifest itself in any macroscopic ther- 
modynamic properties. Application of the same logic to the exact many-body 
eigenvalues would indicate that the many-body localization transition is not fol- 
lowed by any singularities in the static specific heat, etc. On the other hand, at 
this point we can not rule out the possibility that this conclusion is an artifact of 
treating the real parts of the electron self-energies with an insufficient accuracy. 
Most likely scenario, to our opinion, is that the insulating phase behaves like 
a glass (spin or structural) and demonstrates all the glassy properties (Fisher 
and Hertz, 1991; Bouchaud et at, 1998), like absence of ergodicity (even when 
some coupling with phonons is included), effects of aging, etc. Discussion of the 
equilibrium susceptibilities in the latter case becomes quite meaningless. 

The quantitative theory built in BAA paper assumes that the interaction is 
weak. On the other hand, qualitative consideration of the localization of many- 
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body excitations does not rely on this assumption. The important ingredients are 
(i) localization of single-particle excitations, and (ii) Fermi statistics. Consider, as 
an example, Wigncr crystal (Wigner, 1934). It is well known that strong enough 
interaction leads to a spontaneous breaking of the translational symmetry in d- 
dimcnsional clean systems at d > 2. In a clean system Wigner crystallization is 
either a first-order phase transition (d = 3), or a Kosterlitz-Thouless transition 
(d = 2). Even weak disorder destroys both translational and orientational order 
(Larkin, 1970) and pins the crystal. The symmetry of this state is thus not 
different from the symmetry of a liquid, and the thermodynamic phase transition 
is commonly believed to be reduced to a crossover. 

We argue that the many-body localization provides the correct scenario for 
the finite-temperature "melting" transition between the insulating phase, which 
may be called "solid", and the metallic phase, which may be called "liquid". 
Indeed, the conductivity of the pinned Wigner crystal is provided by the motion 
of defects. At low temperatures and in the absence of the external bath, all 
defects are localized by the one-particle Anderson mechanism. Phonon modes 
of the Wigner crystal are localized as well, so the system should behave as a 
many-body insulator. As the temperature is increased, the many-body metal- 
insulator transition occurs. It is not clear at present, whether it occurs before or 
after the crystalline order is destroyed at distances smaller than Larkin's scale. 
Construction of effective theory of such a transition is a problem which deserves 
further investigation. 
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